o 

(N 

U 

Oh 
< 

0\ 



43 

Oh 
i 



CO 

CO 

O 

• *— i 
CO 

43 

Oh 



> 

00 
(N 

O 



X 
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In simulations of a 12.5PW laser (focussed intensity I — 4x 10 23 Wcm -2 ) striking a solid aluminum 
target 10% of the laser energy is converted to gamma-rays. A dense electron-positron plasma is 
generated with a maximum density of 10 26 m~ 3 ; seven orders of magnitude denser than pure e~e + 
plasmas generated with 1PW lasers. When the laser power is increased to 320PW (I = 10 25 Wcm~ 2 ) 
40% of the laser energy is converted to gamma-ray photons and 10% to electron-positron pairs. In 
both cases there is strong feedback between the QED emission processes and the plasma physics; 
the defining feature of the new 'QED-plasma' regime reached in these interactions. 



I. INTRODUCTION 



With the advent of next generation 10PW-100PW 
lasers pQ a new frontier will be reached in high-power 
laser-plasma physics. These lasers will create strong 
enough electromagnetic fields to access strong-field quan- 
tum electrodynamics (QED) processes thought to be re- 
sponsible for cascades of antimatter production in the rel- 
ativistic winds from pulsars and black holes [2]. Strong- 
field QED processes have typically been investigated us- 
ing particle accelerators in experiments arranged such 
that these QED scattering processes can be studied in iso- 
lation [5] (and their cross-sections compared to QED cal- 
culations [1]). This is also true of laser-solid experiments 
where photon and pair production occur in the electric 
fields of the nuclei of high-Z materials far from the laser 
focus 5 . By contrast, the fields in a > 10PW laser's 
focus will cause strong-field QED reactions directly [5J. 
In this case the QED processes strongly modify the ba- 
sic plasma dynamics. Conversely, the rates of the QED 
reactions depend on the electromagnetic fields which are 
determined by the plasma dynamics. As a result of this 
feedback neither the QED nor the plasma physics may 
be considered in isolation, but both must be treated self- 
consistently in the resulting 'QED-plasma' [7]. 

The important strong-field QED emission processes 
are [4]: (1) quantum-corrected synchrotron radiation; 
(2) multiphoton Breit-Wheeler pair production. In (1) 
electrons and positrons in the plasma radiate energetic 
gamma-ray photons when accelerated by the electro- 
magnetic fields of the laser. In process (2) these pho- 
tons interact with the laser fields and generate electron- 
positron pairs. The controlling parameter for these pro- 
cesses is r\ sa j6^I/I s , where / is the laser intensity & 
I s = £qcE 2 J2 = 2 x 10 29 Wcm~ 2 is the laser intensity at 
which the average field in the laser focus is equal to the 
Schwinger field E s = 1.3 x 10 18 Vm -1 ; i.e. the field re- 
quired to break down the vacuum into electron positron- 
pairs [8j. 8 £ [0,2] depends on the interaction geometry. 
As r\ reaches 0.1 (for an optical laser) the electrons in 



the plasma radiate a significant fraction of their energy 
as gamma-ray photons; the plasma enters the 'radiation 
dominated' regime [5]. In this regime the radiation reac- 
tion force [10] must be included in the equation of mo- 
tion of the electrons h, positrons. As laser intensities 
increase from the current maximum of 10 22 Wcm -2 [TT] 
to exceed I ~ 10 23 Wcm~ 2 the ratio I/I s and the 7 fac- 
tor to which the laser accelerates the electrons increase 
in step and r\ — 0.1 is reached. 10PW lasers should be 
able to push well into this regime; if all the energy of a 
10PW laser pulse is focussed into a laser spot of radius 
one micron I = 3 x 10 23 Wcm~ 2 . When rj = 1 the fol- 
lowing quantum corrections to the gamma-ray radiation 
become important: (1) the gamma-ray photon spectrum 
is modified to account for the recoil of the electron as 
it emits [4]; (2) the emitted photon energy becomes a 
significant fraction of the emitting particles energy and 
therefore the emission becomes stochastic [T2]. In addi- 
tion a significant amount of the radiated photons gener- 
ate electron-positron pairs. Therefore when 7/ = 1 the 
'QED dominated' regime is reached [5]. In optical laser- 
plasma interactions this will occur as intensities increase 
to I ~ 10 24 Wcm~ 2 - 10 25 Wcm~ 2 , the limit attainable 
with 30PW-300PW lasers. 

Due to the complexity of the feedback between the 
QED emission processes and plasma physics effects in a 
realistic laser-produced QED-plasma, numerical simula- 
tions of these interactions are essential. The appropriate 
simulation tool is obtained by augmenting a particle-in- 
ccll (PIC) code by including QED emission processes (1) 
& (2) above. A classical model of gamma-ray emission 
and the resulting radiation reaction has previously been 
included in several PIC codes [T3T[l"5] . A classical model 
is only valid in the relatively narrow intensity range defin- 
ing the radiation dominated regime. A quantum treat- 
ment of gamma-ray photon and pair generation, valid 
in the radiation and QED dominated regimes, has only 
recently been coupled to PIC codes [7J UMTS] . 

The resulting QED-PIC codes have been used to self- 
consistently simulate cascades of electron-positron pair 
production, where a critical density pair plasma can 



be generated from a single seed electron, in pulsar at- 
mospheres [17j and in the interaction of two counter- 
propagating 100PW (7 = 3 X 10 24 Wcm~ 2 ) laser pulses 
[15] . QED-PIC simulations have also recently been used 
to show that prolific gamma-ray photon and pair produc- 
tion is possible in fOPW laser-solid interactions [7]. Here 
the interaction of the laser pulse with a dense plasma, 
combined with the fact that the solid reflects the laser 
pulse and so doubles the electric field, compensates for 
the expected QED rate reductions due to the ten times 
lower intensity. In this paper we present a study of the 
interaction of O(10PW) and 0(1 OOPW) laser pulses with 
solid aluminum targets, exploring this promising configu- 
ration for pair production. We will use QED-PIC EPOCH 
[T§] simulations to elucidate the details of the feedback 
between QED emission processes and plasma physics ef- 
fects. In doing so we will attempt to outline a theoretical 
framework for laser-solid interactions in the QED-plasma 
regime, which will be important not only for determin- 
ing the most effective way to produce copious numbers 
of gamma-ray photons and pairs, but also in determin- 
ing the viability of any proposed applications of >10PW 
laser-solid interactions such as ion acceleration to multi- 
GeV energies [201 [2T] and high-harmonic generation [22] . 



II. QED EMISSION MODEL 

In this section we will discuss the model used for the 
QED emission processes its coupling to a PIC code. This 
is dramatically simplified by the fact that in high inten- 
sity laser-plasma interactions the macroscopic laser fields 
are effectively unchanged in the QED interactions [33] . In 
this case the laser's electromagnetic field may be treated 
classically [21] and the QED reactions in the strong-field 
QED framework |25j . Two approximations may be made 
concerning the classical 'macroscopic' laser fields. (1) 
They are quasi-static. The length scale over which pho- 
tons or pairs are formed is a factor of a = eEiAi/2irm e c 2 
(Ei is the electric field of the laser) times smaller than 
the laser wavelength A/ [5B]- For I > 10 23 Wcm~ 2 , a >• 1, 
and the laser's fields may be treated as constant during 
the QED emission processes. (2) The laser's fields are 
much weaker than the Schwinger field. In this case the 
QED reaction rates in the general fields in the plasma are 
the same as those in plane wave fields and depend only on 
the Lorentz-invariant parameters: r\ = (eh/m^c^F^p"] 
and X = {eh 2 2m i e c A )\F lLV k v \ g]. p^ (&**) is the elec- 
tron's (photon's) 4-momentum. For I < 10 25 Wcm~ 2 , 
E/E s < 5 x 1CP 3 and so approximation (2) holds. 

Physically 77 is the field perpendicular to the electron 
motion boosted into its rest frame. This can be seen 
clearly by writing 77 in terms of three- vectors in the ultra- 
relativistic limit: 77 sa (j/E s )\'E± + j3 x cB|; similarly 

X = (fiw 7 /2m e c 2 )|Ej_ + k x cB|. Here Ej_ is the com- 
ponent of the electric field perpendicular to the parti- 
cle's motion (in the direction of k for gamma-ray pho- 
tons and v = /3e/|v| for electrons & positrons). The 



three- vector form for 77 shows the origin of the geometri- 
cal factor 8. In the case of an underdense plasma being 
struck by a single high-intensity laser pulse the electrons 
are rapidly accelerated to ~ c in the direction of propa- 
gation of the laser, the Ej_ and (3 x cB terms in 77 almost 
exactly cancel and so 6 <C 1, f] is small and emission 
of photons and pairs consequently reduced. In the case 
of an ultra-relativistic beam of electrons colliding with 
the laser pulse the terms add and 6 ~ 2, dramatically 
increasing the level of emission. In the case of counter- 
propagating laser pulses 9«1 and the situation is also 
favorable for gamma-ray photon and pair production. A 
similar configuration is found in laser solid interactions, 
with the counter-propagating beam provided by the re- 
flected wave. We will compare the effectiveness of the 
counter-propagating beam and laser-solid configurations 



in producing gamma-rays and pairs in section III A 



A. Emission Rates & Monte-Carlo Model 

For a plane electromagnetic wave the (spin & po- 
larization averaged) rates of photon production by 
electrons and positrons of energy jm e c 2 & pair 
production production by photons of energy Huj 7 
are, respectively: A-y (77) = (v3afc/X c )(ri/^)h(ri) & 
A±(x) = (2na f c/X c ){m e c 2 /h J J 1 )xT ± (x) @]. a f is 
the fine-structure constant, A c is the Compton wave- 



length. h(j]) 



lo d X F(v,x)/x, F(V,X) is the 



quantum-corrected synchrotron spectrum |3]. T±(x) 
0.16K 2 1/3 (2/3 X )/x- 

The quantum (stochastic) nature of the emission |12j 
can be modelled using a Monte-Carlo technique [Tt)ll27| . 
The cumulative probability that a particle emits when 
passing thorough a plasma of optical depth r is P(t) — 
1 — e~ r ; t = J X[rj(t')]dt' and A is the reaction rate. 
To determine the optical depth each particle traverses 
before emitting P is assigned a random value between 
and 1 and the equation for P above is then inverted 
to yield r. For each particle the optical depth evolves 
according to the rate equations above until r is reached 
and emission occurs. The cumulative probability that 
an electron or positron with parameter 77 emits a photon 
with x (i- e - an energy hcu 7 — 2jm e c 2 x/v) is Px( r l>x) — 
[l/h(r])) f Q x dx'F(r], x')/x- When a photon creates a pair 
it is annihilated and its energy shared between the gener- 
ated electron and positron. The cumulative probability 
that the positron takes fraction / of the energy (parame- 

terized by X ) is Pf(f,X) = J f dfp f (f, X ) US]. P x k P f 
are assigned to the emitted photon or positron at random 
in the range [0,1]. The corresponding values of x or / are 
determined by inverting the equations for P x or Pf. Af- 
ter emitting a photon the emitting electron or positron 
recoils, its momentum changing by Ap = —(Hujj/cjp, 
providing the quantum equivalent of the radiation reac- 
tion force PI IM] • Note that in the limit where Ap <§; p 
many photons are emitted for an appreciable change in 



the particles energy and therefore the whole synchrotron 
spectrum is sampled, this is identical to a classical treat- 
ment where the particle instantaneously emits the entire 
synchrotron spectrum, thus the Monte-Carlo algorithm 
agrees with a classical treatment of radiation reaction in 
the classical limit. 



B. QED-PIC 

The basis of the PIC technique 30 is the representa- 
tion of the plasma as macroparticles, each representing 
many real particles such that the number of macroparti- 
cles is amenable to simulation. Particle interactions are 
mediated by: (1) interpolating the charge and current 
densities resulting from the positions and velocities of the 
macroparticles onto a spatial grid; (2) solving Maxwell's 
equations for the E & B fields; (3) interpolating these 
fields onto the particle's positions and pushing the par- 
ticles using the Lorentz force law. The inclusion of the 
QED processes is simplified by: the fact that the macro- 
scopic fields may be treated classically and therefore step 
(2) remains unchanged; the macroscopic fields are quasi- 
static and therefore the QED interactions are point-like, 
occur instantaneously on the timescale of the PIC code 
and are consequently not resolved by the code. Therefore 
we include the QED emission processes as a new step (0) . 

During emission macrophotons and macropairs are cre- 
ated. The pairs are treated in an equivalent way to 
the original electrons in the PIC code. The photons 
are treated as massless, chargeless macroparticles which 
propagate ballistically. The placement of the QED emis- 
sion step at (0) ensures that the feedback defining QED- 
plasmas is simulated self-consistently. Radiation reaction 
exerts a drag force, altering the velocity of the electrons 
and positrons and therefore the current in the plasma; 
pair production acts as a current source. The inclusion of 
the QED processes before the PIC code solves Maxwell's 
equations means that this change in the current is in- 
cluded when the fields are updated. The updated fields 
are then passed back to the QED routines and used to 
calculate emission during the next time-step. 



III. >10PW LASER-SOLID INTERACTIONS 

Two-dimensional QED-PIC EPOCH simulations of 
12.5PW and 320PW laser pulses striking solid alu- 
minum targets at normal incidence have been performed 
which demonstrate the most important aspects of QED- 
plasma physics in both the radiation and QED domi- 
nated regimes. In both cases the laser is linearly p- 
polarizcd and focussed to a spot of radius 1/im on the tar- 
get's surface [i.e. the spatial profile of the laser intensity 
here is ex exp(— y 2 /l/im 2 )]. Temporally the laser power 
P = P for < t < 30fs and P = otherwise. Therefore 
the intensity on-target is 4x 10 23 Wcm~ 2 for P =12.5PW 
and lx 10 25 Wcm- 2 for P =320PW. The target is a fully- 



ionized aluminum foil of thickness 1/im and initial density 
profile p(x, y) =2700kgm~ 3 for < x < L, p(x, y) = 
otherwise. The target is discretized on a spatial grid with 
cell size lOnm and is represented by 1000 macroelectrons 
and 32 macroions per cell (12.5PW case) or 1857 macro- 
electrons and 142 macroions per cell (320PW case). 

Fig. Ilia) shows the results for the 12.5PW laser 
pulse. This shows that prolific gamma-ray (2D blue) 
and positron (red contours) production occur as the laser 
bores into the solid (3D grey). 4.8 x 10 13 gamma-ray pho- 
tons with an average energy of 4.8MeV are produced, cor- 
responding to 10% conversion of laser energy to gamma- 
rays and therefore this interaction is in the radiation 
dominated regime. 10 10 positrons are produced. Despite 
the large number generated the positrons are a minority 
species in the plasma, and therefore the sheath is gen- 
erated by the 'fast' electrons launched into the target 
[31] . In this case the positrons pass through the sheath 
and readily escape the target. A pure electron-positron 
plasma is formed behind the target with a maximum 
positron number density of 10 26 m -3 , 0.1 times the non- 
relativistic critical density for optical lasers. For compar- 
ison the highest positron density outside the target ob- 
tained in 1PW laser plasma experiments is ~ 10 19 m~ 3 . 
It should be noted that similar numbers of pairs are gen- 
erated in each case and that the dramatic increase in 
density is entirely due to the much smaller volume over 
which the pairs are generated in 10PW than in PW laser- 
plasma interactions (~ 1/im 3 compared to ~ 1mm 3 ). 

The average positron energy is 320MeV. This is much 
higher than the average photon energy and approxi- 
mately twice that of the fast-electrons (140MeV). This 
suggests that the positrons are born with relatively low 
energy and are rapidly accelerated by the laser to an 
energy equal to that of the fast-electrons and further 
accelerated by the sheath fields on leaving the target. 
The sheath field acts to confine the fast-electrons inside 
the target and so accelerates positrons [5], doing work 
<frm e c 2 approximately equal to the fast electron energy 
|32) . In this case we expect the average Lorentz fac- 
tor of the positrons to be (7) w a^ B + $ ~ ^ a HB = 
2eEj? l B \i HB /2irm e c 2 « 300MeV, which is consistent 
with the simulations. In total 0.01% of the laser energy is 
converted to positron energy and so their relative effect 
on the plasma dynamics is small. 

Figure [lib) shows simulation results for a 320PW laser 
striking a solid aluminum target. 10 16 gamma-ray pho- 
tons and 10 13 positrons are produced with average en- 
ergies of 92MeV & 2.2GeV respectively. The maximum 
positron density is 1.8 x 10 30 m -3 , an increase of four or- 
ders of magnitude for only a factor of 25 increase in laser 
intensity. 40% of the energy is converted to gamma-rays 
and 10% to electron-positron pairs. Therefore at this 
extreme laser intensity both gamma-ray photon and pair 
production are crucial to the plasma dynamics and there- 
fore the interaction is in the QED dominated regime. 



(a) 




Laser 




Figure 1. (Color online) (a) Gamma-ray photon (2D blue) & positron (red contours) production in the interaction of a 12.5PW 
laser pulse with a solid aluminum target (3D grey), (b) Equivalent plot for a 320PW laser pulse striking solid aluminum. 



The Effect of Plasma Physics Processes on the 
QED Rates 



The key feature of a QED-plasma is the feedback be- 
tween QED and classical plasma physics processes. In 
this section we discuss the effect plasma physics pro- 
cesses have on the rates of the QED reactions in laser 
solid interactions. This is best illustrated by compar- 
ing pair production in the interaction of a laser of inten- 
sity / with a solid target to the alternative configuration 
consisting of two counter-propagating laser pulses of in- 
tensity 1/2 interacting with a low-density gas [SI HE] as 
in the laser-gas case complicated plasma effects are less 
important. The laser-solid configuration has the clear 
advantages that the peak electric field is double that of 
laser-gas case due to reflection and that the pulse inter- 
acts with a dense plasma so that many pairs and photons 
may be produced even when the rates of reaction are low. 

A parameter scan of the effect of increasing laser in- 
tensity on the number of pairs produced by each config- 
uration is conducted using one-dimensional EPOCH sim- 
ulations. The targets considered are: solid aluminum 
(density 2700kgm -3 ) and a hydrogen gas-jet (density 
0.02kgm~ 3 ). The solid targets are semi-infinite to avoid 
complications caused by the laser pulse breaking through 
the target. Figure [2ja) shows the number of positrons 
produced by each configuration. Due to the advan- 
tages previously mentioned the laser-solid case produces 
more positrons for / < 8 x 10 23 Wcm~ 2 . For / > 
8 x 10 23 Wcm~ 2 the gas-jet configuration continues to 
behave as expected, with increasing intensity leading to 
increased pair production and when r\ ~ 1 a large frac- 
tion of the pairs generated go on to produce additional 
pairs, the reaction runs away and a cascade of antimatter 



production ensues. This is in good agreement with the 
results of Nerush et al p~8] . In contrast pair production 
in the laser-solid case peaks when I = 8 x 10 23 Wcm -2 
and then decreases for further increases in laser intensity. 

The difference between the laser-solid and laser-gas 
configurations is more marked when considering the rate 
of pair production, shown in Figure [2jb) . The rate is 
substantially lower for the solid than the gas target at 
all intensities. Several plasma effects have been pro- 
posed as being responsible for this reduction, namely: 
relativistic hole boring, the skin effect & relativistic 
transparency [71 113|. Qualitatively the reduction when 
/ < 8 x 10 23 Wcm~ 2 is due to hole-boring & the skin ef- 
fect. When the pulse strikes the solid surface its radiation 
pressure accelerates the surface to speed v^g = Phbc in 
a process known as hole-boring. The laser is reflected 
in the rest frame of the surface, in which its intensity 
is reduced by a factor of (1 — (3hb)/(1 + Phb), where 
0hb = vS/(l + vS) and H = I/pc 3 is the pistoning pa- 
rameter for a laser of intensity / and a target of density 
P [2J]- From this formula one can see that in the 12.5PW 
interaction vhb ~ 0.2c and the intensity in the rest frame 
of the surface is 0.7 times that in the lab frame. Th elec- 
tric field of the laser is evanescent inside the overdense 
solid and is reduced to E s ° l B = 2{n c /n eHB ) 1 l 2 E™ B x (the 
skin effect). n c = jm e eouj' 2 / e 2 is the relativistically cor- 
rected critical density for the plasma (7 is the Lorentz 
factor to which the electrons are accelerated by the laser 
pulse) ; E BB X is the peak laser electric field and n e HB the 
electron number density both in the hole-boring frame. 
The reduction in the rate of pair production is consis- 
tent withe the reduction in the field in the solid target 
to Ejj B [7J . The Lorentz factor reached by the electrons 
in the solid is oc yl, therefore if the laser intensity is in- 
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Figure 2. (Color online) (a) Number of positrons generated 
in the interaction of a laser pulse of intensity I with solid and 
gas targets, (b) The rate of positron production in each of 
these cases. 



creased eventually n c > n e and the solid becomes under- 
dense and therefore transparent [33] - This occurs when 
I > 8x 10 23 Wcm~ 2 . In this case the electrons are pushed 
forwards at c, the situation is similar to an underdense 
plasma illuminated by a single laser pulse and emission is 
drastically curtailed. Therefore gamma-ray and pair pro- 
duction are maximized in laser-solid interactions when 
the solid is marginally overdense. When this is the case 
the ratio n e /n c is minimized and this overcomes the effect 
of the increased reduction of Ehb due to the increase in 
vhb (resulting from the decrease in p and corresponding 
increase in S). This was shown with simple quantitative 
estimates in [7] . However some aspects of the model used 
there are inconsistent with the work in [14 and further 
work is required. 

The trend of decreasing levels of pair production with 
increasing intensity does not continue. Prolific pair pro- 
duction was seen in the 320PW laser-solid interaction, 
despite the fact that at the intensity reached in this 
interaction (10 25 Wcm~ 2 ) the solid target is relativisti- 
cally underdense. This could be due to a new QED- 
mediated laser absorption mechanism which operates in 
underdense plasmas proposed by Brady et al 34]. How- 
ever care should be taken when comparing to the results 



in this paper, which discussed the radiation dominated 
regime. In the QED dominated regime reached in the 
320PW laser-solid interaction pair production is expected 
to lead to the generation of critical density pair plasmas 
over the duration of the laser pulse [18| which will clearly 
influence the plasma physics and as a result the macro- 
scopic fields and therefore the QED rates. For example, 
it has recently been shown that in this regime a dense 
pair-plasma can form in front of the solid surface (35) . 
This pair plasma can absorb the laser. In this case we 
expect the QED rates in the solid to be reduced. Much 
more work is required to fully understand the influence of 
plasma physics processes on the emission rates in laser- 
solid interactions in the QED dominated regime. 



B. The Effect of QED Emission on the Plasma 
Physics 

The substantial amount of energy converted to 
gamma-ray photons and pairs profoundly alters the laser 
energy absorbed by the electrons and ions in a QED- 
plasma and so the plasma processes which are driven 
by this energy. First we consider the interaction of the 
12.5PW laser pulse with the solid aluminum target. In 
this case 10% of the laser energy is converted to gamma- 
ray photons. The effect that this has on the energy spec- 
tra of the electrons and ions is shown in Figures p[a) & 
(3tb). The average energy of the fast electrons drops from 
150MeV to 140MeV. It is clear that the most energetic 
electrons are most affected as they emit photons most 
strongly. Gamma-ray emission and the resulting radia- 
tion reaction causes a substantial difference in the ion 
spectrum, which develops two peaks. Preliminary qual- 
itative discussions of the modification of some plasma 
processes caused by this change in the energy spectra are 
given in Refs. |14) . 

Next consider the I = 1 x 10 25 Wcm -2 laser-solid inter- 
action. This interaction is in the QED dominated regime 
and a significant fraction of the laser energy is converted 
to both gamma-ray photons and pairs (40% & 10% re- 
spectively). Figure |3J^c) shows the effect QED processes 
have on the electron energy. The average fast electron 
energy is reduced from 4.5GeV to 2.0GeV. As before, 
the high energy tail is preferentially damped. However, 
in this case pair production is a significant source of elec- 
trons, which are generated at moderate energies and sig- 
nificantly enhance the spectrum here. Figure p[d) shows 
that QED effects do not significantly alter the ion spec- 
trum. In this underdense case the ions gain energy by 
coupling to electrostatic fields generated by the electrons 
as they are pushed forwards, as we have seen emission 
is not strong for electrons undergoing such motion. We 
expect the dramatic modification of the electron and ion 
spectra caused by QED emission in the QED dominated 
regime to strongly effect the plasma processes; however, 
very little work has been done to elucidate this. 
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Figure 3. (Color online) Energy spectra of the plasma components with and without the inclusion of QED effects in both the 
12.5PW and 320PW laser-solid simulations, (a) & (c) show electron spectra for 12.5PW and 320PW respectively, (b) & (d) 
ion spectra for 12.5PW & 320PW. 



IV. CONCLUSIONS 



Next generation high-power lasers, operating at in- 
tensities > 10 23 WcnU 2 , will generate a qualitatively 
new plasma state on interacting with matter. These 
QED-plasmas are defined by feedback between QED 
emission processes and classical plasma physics effects. 
We have described how the important QED processes: 
synchrotron-like gamma-ray photon emission & multi- 
photon Breit- Wheeler pair production; can been included 
in a PIC code and used the resulting QED-PIC code to 
simulate this feedback in >10PW laser-solid interactions 
self-consistently. We have shown that the rates of re- 
action in the simulations can only be explained when 
plasma effects are included and that the QED modifies 
these plasma effects by strongly altering the electron and 
ion energy spectra. This alteration of the energy budget 
of laser-solid interactions may be important for proposed 
applications of 10PW lasers, such as ion acceleration or 



harmonic generation. Simulation of a 12.5PW laser pulse 
striking a solid aluminum target demonstrates the con- 
version of a significant fraction (10%) of the laser en- 
ergy to gamma-ray photons. In addition a pure electron- 
positron plasma is generated in the simulation with den- 
sity seven orders of magnitude higher than currently 
achievable in laser-matter interactions. Simulations of 
a 320PW laser-solid aluminum target interaction demon- 
strate that in this case we expect not only efficient (40%) 
conversion of laser energy to gamma-rays, but also (10%) 
to pairs. This prolific production of gamma-ray photons 
and pairs may find application as an efficient and bright 
sources of these particles. 
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